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Abstract 

Given a rational algebraic curve denned by means of a rational parametrization, 
we address here the problem of detecting whether the curve exhibits some kind of 
symmetry (central, mirror, rotation), and of computing the elements of the symme- 
try in the affirmative case. We provide effective methods for solving these questions 
without any conversion to implicit form. The underlying idea is the existing rela- 
tionship between two proper parametrizations of a same curve, which in turn leads 
to algorithms where only resultants and univariate gcd's are involved. These meth- 
ods have been implemented and tested in the computer algebra system Maple 15; 
evidence of their applicability, as well as a detailed theoretical analysis, are given. 



1 Introduction 

The problem of detecting the symmetries of an algebraic curve has been ex- 
tensively studied mainly because of its applications in Pattern Recognition. A 
common problem in this area is how to choose, from a database of algebraic 
curves, the one which best suits a given object, also represented by means 
of an algebraic equation (see for example [6], [10], [15], [18], [19], [20]). For 
this purpose, first one must place the shape to be identified in a "canonical 
position", so that the comparison can be carried out. Hence, as a main step, 
the symmetries of the studied curve must be computed. Among others, the 
computation of symmetries has been addressed in [6], using splines, in [2], 
[3], [22], by means of differential invariants, in [8], [9], [18], using a complex 
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representation of the implicit equation of the curve, or in [6], [16], [17], [21], 
using moments. 



In this paper we address the problem of computing the symmetries of a 
real plane curve C defined by means of a rational parametrization <p(t) = 



a consequence, our input is exact. More precisely, we consider three types of 
symmetry: central symmetry (i.e. symmetry with respect to a point, which 
is the center of gravity of the curve), mirror symmetry (i.e. symmetry with 
respect to an axis), and rotation symmetry (which means that the curve is in- 
variant under a non-trivial rotation around a point). Notice that while central 
symmetry is a particular case of rotation symmetry (namely, with 6 — 7r), it is 
interesting enough as to be addressed separately. In this context, we provide a 
novel approach which allows to detect whether C exhibits some kind of sym- 
metry or not, and to compute the elements of the symmetry in the affirmative 
case. Our method provides a deterministic yes/no answer to the question on 
the existence of a certain kind of symmetry. As for the computation of the el- 
ements of the symmetry, since our technique requires to numerically compute 
the real roots of certain polynomials, in general the output is also numerical. 
However, a numerical test can be applied to evaluate the exactness of the 
computed element (see Section 7). 

So, our paper differs in two ways from other papers on the subject. On the 
one hand, while most papers assume that the curve to be analyzed is given by 
means of its implicit equation, we work with the curve directly in parametric 
form. In fact, our method does not require to compute the implicit equation 
of the curve. This is an advantage whenever the input is defined in parametric 
form, which is often the case, for example, in Computer Aided Geometric De- 
sign (CAGD). On the other hand, almost all the papers address the question 
directly in a numerical way, and as a consequence do not provide in general de- 
terministic answers to the question of whether the curve is symmetric or not. 
This is a reasonable strategy when one works with fuzzy objects (i.e. images), 
which is often the case in Pattern Recognition. However, it is not completely 
satisfactory when one deals with an exact object, and therefore is expecting 
also an exact answer. Still, some papers in the literature have considered some 
of these two questions. For the first one, parametric curves have been consid- 
ered in [6], where the curve is assumed to be a spline (i.e. a union of pieces of 
polynomial parametrizations), and in the papers on differential invariants ([2], 
[3] and others), where the input is considered to be a parametrization without 
any restriction on its functional form. Nevertheless, in [6] the algorithm which 
is provided is essentially numerical. In the case of the papers on differential 
invariants, even though in principle the idea is applicable to exact inputs, the 
subsequent computations turn too hard unless the input has a very low de- 
gree. In fact, in these papers the goal is to translate the idea into numerical 




where Pi(t),qi(t) are polynomials with integer coefficients; as 
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methods, which do work efficiently. As for the second question, in [8], [9] we 
find algorithms on implicit curves (i.e. the implicit equation must be known) 
with a good performance over exact inputs; however, here the method requires 
the leading coefficients of the implicit equation of the curve to be different, 
which is a limitation in certain cases. 

The main ingredient behind our method is the existing relationship between 
two proper (i.e. injective up to finitely many parameter values) parametriza- 
tions of a same curve. The rough, general idea is the following: if the curve 
C has a certain symmetry, then by applying this symmetry we can obtain 
another parametrization <f>(t) from the original 4>(t); whenever (f>(t) is proper, 
then <j)(t) will be also proper, and both parametrizations will be related by 
means of a certain transformation whose general form is well-known (it is a 
Mdebius transformation) , depending on 4 real parameters. This argument can 
be reverted as to characterize the existence of the symmetry. Hence, we detect 
symmetry iff we can find real values for these parameters. For each kind of 
symmetry, we can show that these 4 parameters can be expressed as rational 
functions of just one of them. Thus, in the end the problem boils down to 
computing greatest common divisors of univariate polynomials, and checking 
whether certain univariate polynomials have some real root, or not. In this 
sense, we have implemented our algorithms in the computer algebra system 
Maple 15, and we have tested them over many examples. We report on this 
in the last section of the paper. Additionally, we also present an improve- 
ment of our method for curves which admit a parametrization of the type 



r > or s > 0. This family is important from the point of view of applications 
because it contains most of the trigonometric curves, i.e. curves which are 
parametrized by truncated Fourier Series (see for example [1], [5], [12], [13]). 

The structure of the paper is the following. In Section 2 we present some pre- 
liminary results to be used throughout the paper. The detection of central, 
mirror and rotation symmetry is addressed in Section 3, Section 4 and Section 
5, respectively. The special type of parametrizations mentioned above is ad- 
dressed in Section 6. Finally, implementation issues are considered in Section 
7. An appendix contains the (long) proof of a result which is needed in Section 
4. 



2 Preliminary Results 

Along the paper we consider a plane algebraic curve C defined by means of a 




where p(t),q(t) are polynomials, r,s 



€ N and either 



rational parametrization <p(t) = (x(t),y(t)) where x(t) 



y(t) 
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with pi, qi G R[t] for i — 1,2, and gcd(pj, = 1 for 2 = 1, 2. Furthermore, we 
will assume that 0(0 is proper, i.e. that it is injective for almost all (complex) 
values of t. Algorithms for checking properness can be found for example in 
[14]. Also, from the algorithm in page 193 of [14] it follows that every rational 
curve can be properly reparametrized without extending the ground field. Now 
the following result will be crucial for us. 

Theorem 1 Let 1 (£),0 2 (O G R[t] be two proper rational parametrizations 

of a same curve. Then there exists a unique function (pit) = — — , with 

7* + 

a5 — /?7 7^ 0, a, (3, 7, 5 G R, fulfilling that fait) = 0i (</?(£)). 

Proof. The existence of a, (5, 7, 5 G C fulfilling 02 (0 = 0i(y ; >(O) i s guaranteed 
by Lemma 4.17 in [14]. So, we just have to prove that a, (3, 7, 5 can be assumed 
to be real. For this purpose, observe first that tp(t), if it exists, is unique. 
Indeed, if ipi(t), y?2(0 both fulfill the conditions in the statement, we get 
that 02(0 = 0i(v?i(O) = 0i(<£2(O)- Since 0i(O is proper, we deduce that 
(pi(t) = (p2{t) for almost all complex values oft. Furthermore, since (pi(t), v?2(0 
are analytic, ipi(t) = v?2(0 follows from the Identity Theorem (see e.g. page 81 
in [7]) Thus, in order to prove the statement, we just need to find one ip(t) with 
real coefficients. Now from page 97 in [14] one may see that (p(t) = (f)^ 1 o 2 
fulfills the conditions of the statement. Also, it has real coefficients because 
0i (O>02(O are parametrizations with real coefficients, and the computation 
of the inverse 07 1 does not extend the ground field (see page 107 in [14]). n 

The functions (p(t) = — - are called Mdebius transformations. In the rest 

7* + 

of the paper we will maintain the notation ip(t) for this type of functions, and 
we will assume, according to the above result, that a, (3, 7, 5 G R. Also, in the 
sequel we will often use the complex notation z(t) = x(t) + iy(t) for the points 
of C, therefore seen as elements of C. Furthermore, in the rest of the paper we 
will denote as || ■ || the Euclidean norm over C. Additionally, along the paper 
we will speak of the degree of the parametrization 0(0 to denote the maximum 
power in the numerators and denominators of its components. 

In the sequel, we will assume that C is neither a line nor a circle, which can 
be considered as trivial cases from the point of view of symmetry analysis. 



3 Central Symmetry 



In [8], [9] it is proven that the center of symmetry of an algebraic curve, if it 
exists, is unique. Furthermore, the curve C has central symmetry with respect 
to a point z (in complex notation) iff z — (z(t) — z ) = 2z — z(t) is also 
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a point of C for every value of t (notice that 2z — z(t) corresponds to the 
symmetric of a generic point z(t) of C with respect to z ). This happens iff 
z(t) = 2z — z(t) is also a parametrization of C. Furthermore, if z(t) is proper 
(which holds by hypothesis) then z(t) is also proper. So, the following result 
follows from Theorem 1. 

Theorem 2 The curve C is symmetric with respect to a point zq iff there exist 
a, (3, 7, 5 e R, a5 — f3j ^ 0, such that 

z(<p(t)) = 2z - z(t) (1) 

Hence, detecting central symmetry is equivalent to check whether there exist 
(p(t), zq in the conditions of Theorem 2. In turn, these conditions lead to the 
analysis of a polynomial system (with six unknowns: the four parameters in- 
volved in tp{t) and the two coordinates of z ). In general this system is too 
complicated. So, let us find extra conditions on the parameters of <f(t). We 
start with the following lemma. Here, we denote <p 2 (t) = (ip o (p)(t). 

Lemma 3 Assume that C is symmetric with respect to z . Then the function 
(p(t) satisfies that (p 2 (t) = t. 

Proof. From z(ip(t)) = 2z — z(t), we get that z((p 2 (t)) = 2z — z((p(t)) = 2z — 
(2z — z(t)); so, z(ip 2 (t)) = z(t). Now since z(t) is a proper parametrization, 
we get that p 2 {t) = t for almost all (complex) values of t, i.e. <f 2 (t) — t = for 
infinitely many (complex) values. Since <p 2 (t) — t is analytic, from the Identity 
Theorem (see page 81 in [7]) it follows that it must be identically zero. □ 

By explicitly computing <p 2 (t), we deduce the following corollary from the 
above lemma. 

Corollary 4 IfC is symmetric with respect to z , then (a+S)j3 = 0, = 
0,a 2 = 5 2 . 

Let us distinguish now two cases for <p(t), namely 5 = and 5^0, which 
are analyzed in the following subsections. Also, for technical reasons in the 
sequel we will assume that z(t) (and therefore its derivatives) is well defined 
for t — 0; notice that this can always achieved by applying a linear change of 
parameter (which does not affect properness of the curve). 

3.1 Case 5 = 

From Corollary 4, a 2 = 5 2 . So, a = and ip(t) = k/t, k G M. Now if C has 
central symmetry and ip(t) has this form, by differentiating the equality (1) 
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we get 



z'(k/t) ■ {-k)/t 2 = -z'(t) (2) 
Lemma 5 The equation (2) cannot be an identity \/k. 

Proof. Substituting k = in the above equation, we deduce that z'(t) is 
identically 0, which implies that z{t) constant. So, C would not be a curve, 
which is absurd. □ 



Let be the gcd of the numerator in (2) (considered as a polynomial in 
t). From the above lemma, we have that £(k) is not identically 0. Hence, we 
deduce the following result. 

Theorem 6 If C has central symmetry and <p(t) = k/t, then £(k) = 0. Con- 
versely, if k G R is a real root of and the denominator in (1) does not 
vanish, then C has central symmetry with <p(t) = k /t, and the symmetry 
center can be computed from (1). 



3.2 Case 8^0 



Since a 2 = 8 2 , we can distinguish the subcases a — 8, and a = —8, respec- 
tively. In the first case, since we are assuming 5 / we have a + 8 ^ 0, and 
from Corollary 4 we get (3 = 7 = 0. So, <p(t) = t. However, in this situation 
from Theorem 2 we get that z(t) is a constant, i.e. C is not a curve. So, the 
only possibility is a = —8 7^ 0. Now since 8 7^ 0, by dividing if necessary the 
numerator and denominator of ip(t) by 8, we can assume that 8 = 1; thus, 
a = —1. Hence, if C has central symmetry then by differentiating (1) and 
evaluating at t — 0, we get 

z'((3). (1 + P 1 ) = z'(0) (3) 



We distinguish two cases, depending on whether f3 = or f3 7^ 0. If (3 = 0, 
then ip(t) = — - . In this situation, we consider the equation 

z'( V (t)) V '(t) + z'(t) = (4) 



which results from differentiating (1) and substituting ip(t) = . One 

7* + 1 

may argue as in Lemma 5 to derive that this equation cannot be an identity. 
So, the gcd of the numerator 2/(7) of the above expression (considered as a 
polynomial in t) cannot be identically 0. If this polynomial has some real 
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root, by substituting t — in (1) we get that C has central symmetry, with 
symmetry center z(0). Hence the following result holds. 

Theorem 7 If C has central symmetry and (pit) = — — with 5 ^ 0, (3 = 0, 

yt + 5 

then (i) a = — 1, 5 = 1; (ii) v( r ))=0. Conversely, if 70 € M. is a real root of 
1^(7), snc/i #ierf no denominator in (1) vanishes identically, thenC has central 

symmetry with ip(t) = , and the symmetry center is z(0). 

If P 7^ 0, from (3) we get 



7 



V(0) 



1 

i9 



The right hand-side of the above equality is a complex number, but from 
Theorem 1 we know that 7 can be taken real. So, denoting 7 = /(/?) + ig(/3) 
we have that 7 = /(/?) and o(/3) = 0, where f{0),g{0) are rational functions 
of/3. 

Lemma 8 Under our assumptions, g(f3) is not identically 0. 

—x'(0)y'(P) + y'(0)x'((3) 

Proof. A direct computation shows that g(j3) = 77777, ^777 • So, 

x' 2 (f3) + y' 2 (/3) 

g(/3) is identically iff x'(0)y'(/3) = y'(0)x'(/3) for all f3 e M. However, in such 
case C is a line, which is excluded by hypothesis. □ 

Finally, consider the equation (4) with ip(t) = JjJj^TJ^i ' us denote by 

the gcd of the coefficients in t of the numerator of this expression. Also, 
let rj(/3) = gcd(g(/3), £(/?)). Then the following theorem holds. 

Theorem 9 If C has central symmetry and ip(t) = — - with 5^0, then 

'jt + d 

(i) a = —1, 5 = 1, 7 = f((3); (ii) i](f3) = 0. Conversely, i/ ft 6 K is a real 
root ofrj((3), such that f(/3) is well defined and no denominator in (1), (3), 

vanishes identically, then C has central symmetry with ip(t) = jy-^ yr~~7r , an d 

the symmetry center can be computed from (1). 



3.3 Full Algorithm 



The following algorithm follows from the preceding subsections. 
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Algorithm Central Symmetry: Given a curve C by means of a proper 
parametrization z(t) (in complex form), well-defined at t — 0, the algorithm 
checks if it has central symmetry, and computes the symmetry center 
in the affirmative case. 

(1) (tp(t) = k/t) Check if has any root ko G K in the conditions of Theorem 
6. In the affirmative case, compute the symmetry center zo, and return 
"The symmetry center is Zq" . 

(2) (^p(t) = + ^ ^ °' ^ = °) Check if v (l) nas an y root 7o ^ M in the 
conditions of Theorem 7. In the affirmative case, return "The symmetry 
center is z(0)" . 

(3) (<p(t) = ~1 , 5 ^ 0, /3 ^ j Check if rj(/3) has any root f3 e M in the 

conditions of Theorem 9. In the affirmative case, compute the symmetry 
center z Q , and return "The symmetry center is Zq" . 

(4) If (1), (2), (3) have not succeeded, then return "The curve has not central 
symmetry" . 

Lemma 5 and Lemma 8 guarantee that the algorithm terminates. Its correct- 
ness follows from the preceding results in this section. 



4 Mirror Symmetry 



Let L be a line passing through a point zq (in complex notation) and forming 
an angle 9 with the x-axis. Then C exhibits mirror symmetry with respect to 
C iff the curve C, parametrized by (z(t) — zo)e ld , obtained by applying over C 
a translation of vector — Zq followed by a rotation around the origin of angle 
—9, is symmetric with respect to the This condition is equivalent to 

(z(t) — zo) ■ e~ ie being also a point of C for every value of t. So, from Theorem 
1 the following result follows. 

Theorem 10 The curve C is symmetric with respect to a line C passing 
through a point zq (in complex notation) and forming an angle 9 with the 
x-axis iff there exist a, (3, 7, 5 G R, a5 — f3j 7^ 0, such that 

[z(<p(t))-z }-e 2ie = zJtj-zo (5) 

Lemma 11 Assume that C is symmetric with respect to z . Then the function 
ip(t) satisfies that (p 2 (t) = t. 

Proof. Substituting t by tp(t) in (5), we get that [z(^ 2 {t)) — z o ]-e 210 = z(tp(t)) — 
Zq. Substituting here z(ip(t)) also from (5), we get z((p 2 (t)) = z(t). Now one 
argues as in Lemma 3. □ 
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So, Corollary 4 holds also in this case, and we can consider the same discussion 
as in the section before. For this purpose, in the rest of the section we will 
assume that z(0) (and therefore all the derivatives z^(0), for p e N) are 
well-defined. This can be achieved by almost all linear changes of parameter. 
Additionally, w.l.o.g. we will suppose that z'(0) ^ 0. 



4.1 Case 5 = 

Here ip(t) = k/t. Now differentiating (5) w.r.t. t, we get 

z'(k/t) ■ (-k) -e m = V{t) -t 2 (6) 

In particular, this implies that 



z'(k/t)-(-k) K) 



Differentiating (6) again, 

z"(k/t) ■ k 2 ■ e 2W = z^j ■ t A + At) ■ 2t 3 (8) 

By dividing (6) and (8), and clearing denominators, we get that 



- z\k/t) ■ [z"{t) ■ t 2 + z'(t) ■ 2t] = z"(k/t) ■ k ■ z'(t) (9) 



In addition to this, taking modules in (7) it holds that 

7(t) -t 2 



z'(k/t) ■ (-k) 



(10) 



Lemma 12 The equations (9) and (10) cannot be identities at the same time 
Vk. 

Proof. Assume that both expressions are simultaneously identities. In that 
case, (9) is a differential equation, and by integrating it we get t 2 z'(t) = 



C • z'(k/t). By taking (10) into account, we deduce that C = ±k. So, t 2 z'(t) = 
±k- z'(k/t). Since this last equality holds for all t, k, making k — (recall that 
z'(0) is well defined) we get that z'(t) = 0; so, z(t) is constant and therefore 
C is not a curve, which is a contradiction. n 

Now let us represent by £i(fc) the gcd of the coefficients in t of the polynomial 
in (9), and let ^(k) be the gcd of the numerator of the polynomial obtained 
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when clearing denominators in (10). Also, let 



d_ ( z'jt) ■ t 2 
Jt \z'(k/t) ■ (-k) 



) 



= 



(11) 



and let £ 3 (k) be the gcd of the coefficients in t of the numerator of the 
above expression. Notice that the condition £ 3 (k) = expresses the prop- 
erty that the right-hand side of (7) is a constant complex number; this con- 
dition, together with the condition ^(k) = 0, ensures that the right-hand 
side of (7) is a complex number of modulus 1, and therefore that (7) holds. 
Now if C exhibits mirror symmetry with (p(i) = k/t, then k must be a com- 
mon root of (,i(k) , t;2(k) , £,s(k) . On the other hand, by Lemma 12, i](k) = 
gcd(^i(k),^ 2 (k),^(k)) cannot be the zero polynomial. 

Theorem 13 // C has mirror symmetry and ip(t) = k/t, then i](k) = 0. 
Conversely, if k e K fulfills i](ko) = 0, and the denominators in (5), (7), 
(9), do not vanish, then C has mirror symmetry and the symmetry axis can 
be computed from (5). 

The equation of the symmetry axis is obtained after performing the substi- 
tution of (p(t) and e 2l ° (in terms of k,t) on (5). Indeed, this way we reach a 
complex expression 

Az + Bz^ + C = 

with A, B, C being complex numbers. Regarding z as a complex variable, and 
substituting z = x + iy, the equation of the real line is the gcd of the real and 
imaginary parts. 

4.2 Case 5^0 

Arguing as in Section 3, we distinguish the subcases a — 5, and a = —5. In 
the first case, we conclude that (p(t) = t. However, let us see that in this case C 
exhibits mirror symmetry only if it is a line. Indeed, if tp(t) = t and C has this 
type of symmetry, then from Theorem 10 we have that (z(t) — Zo)e ld = (z(t) — 
z^)e~ ie . So, (z(t) - z )e ie = (z(t) - z )e w , and hence (z(t) - z )e ie G R. One 
may check then that the imaginary part of (z(t) — z )e ld is (x(t) — x )sin(#) + 
(y(t) — y )cos(6), and this expression must be identically 0. So, we get that 



y(t) - yo 

x(t) - x 



tg(0) = constant 



which means that C is a line. 
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Hence, we only need to consider the case a = —5. In this case, since 5 7^ 0, we 
can assume 5 = 1, and so a = — 1. Now differentiating (5) with respect to t, 
we get 



z'(<p(t)) ■ <p'(t) ■ e™ = z'{t) 



(12) 



where ^ = Jt^TW 

t = we get that 



, and A = a5 — f3j 



z'{/3) ■ A • e 2ie = z'(Q) 



f3j. So, evaluating (12) at 



(13) 



Also, from here 



J2\6 



z'(/3) ■ A 



(14) 



Differentiating again (12) and evaluating at t = 0, we get 

z"(f3) ■ A 2 • e m = z 17 ^) + 2 7 7(0) 



(15) 



By combining the above two equations, we deduce that 

-1 
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z\(5) ■ £22 + z"{P) 
z'(0) 



f3z»(f3) + 2z'(P) 



(16) 



The above expression does not make sense when tz"(t) + 2z'(t) is identically 
0. However, this corresponds to a differential equation that can be integrated, 

whose solution is z(t) = C± H . And it is easy to see that this corresponds 

to a line, which is a case excluded by hypothesis. Now let us write 7 = 7(/3) = 
/(/3) + i£i(/3). Since 7 can be taken real, we have that 7 = f(/3), and = 0. 
Thus A = — 1 — f3f(/3); by taking modules in (14) we get that 



\ zW - « 







(17) 



Now let us represent by ^(/S) the numerator in (17). Then the following result 
holds. The (quite long) proof of this result is provided in Appendix I, so as 
not to stop the flow of the paper. 

Lemma 14 Under our hypotheses, if C exhibits mirror symmetry, then the 
polynomials and ^(P) cannot be identically at the same time. 
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Additionally, let us consider the equation 



d ( z'(t) 

It \*{<p{t)) ■ <p'{t) 



) 







(18) 



which is obtained by solving (12) for e 2ld , and differentiating the resulting ex- 
pression. This condition, together with (17), guarantees that (12) holds. Let 
£ 3 (/3) be the gcd of the numerator of (18), and let rj(/3) = gcd(£i(/3), ^(P), £3 (/?)), 
which, from Lemma 14, is not identically 0. Then we have the following the- 
orem on the existence of mirror symmetry. 

Theorem 15 If C has mirror symmetry and (p(t) = - ^ — with 5^0, then 

it holds that: (i) a = —1, 8 — 1, 7 = f{0); (H) = 0- Conversely, if: (i) 
f3 E R fulfills fj((3o) = 0; (ii) no denominator in (5), (14), (15), (16), (18) 

vanishes identically, then C has mirror symmetry with ip(t) = ^ and 

the symmetry axis can be computed from (5). 



4.3 Full Algorithm 

The following algorithm follows from the preceding subsections. 

Algorithm Mirror Symmetry: Given a curve C by means of a proper 
parametrization z(t) = x(t) + iy(t) (in complex form), well-defined at t — 0, 
such that z'(0) 7^ 0, the algorithm checks if it has mirror symmetry, and 
computes the symmetry axes in the affirmative case. 

(1) (ip(t) = k/t) Check if rj(k) has any real roots k\, . . . , k r e K in the condi- 
tions of Theorem 13. In the affirmative case, let Ci, . . . , C r be the symmetry 
axes they correspond to. 



the conditions of Theorem 15. In the affirmative case, let Ci,...,C s be the 
symmetry axes they correspond to,. 
(3) If (1) and (2) have not succeeded, then return "The curve has not central 
symmetry". Otherwise, return the list of £,'s and £/s. 

Lemma 12 and Lemma 14 guarantee that this algorithm terminates. The cor- 
rectness follows from the results in this section. 




Check if r)(/3) has any real roots j3i,...,/3 s in 



12 



5 Rotation Symmetry 

In [8], [9] it is proven that the center of rotation of an algebraic curve, if 
it exists, is unique. Furthermore, the curve C has rotational symmetry iff 
there exists a point Zq and an angle 9, such that for every value of t, z(t) = 
zo + e l6 * ' (z(t) — zo) is also a point of C (notice that this expression describes 
a rotation of z(t) around zq). Hence, from Theorem 1 we get the following 
result. 

Theorem 16 The curve C has rotational symmetry with center zq and angle 
9 iff there exist a, (3,^,5 G R, a5 — f3j ^ 0, such that 

z(<p(t)) = z + e w ■ (z(t) - z ) (19) 

In [8] it is proven that 9 has the form 9 — — , where n G N. Moreover, from 
Bezout's Theorem one may see that n < 2d, where d is the degree of C. Now 
in the sequel let us assume that C has rotational symmetry, and let us find 
conditions on the parameters of (pit). We will suppose that z(0) is well-defined, 
and z'(0) ^ 0. As in the other sections, we distinguish the cases 5 = and 
5 7^ 0, respectively. 



5.1 Case 5 = 

We consider first the following lemma, which can be proven in a similar way 
to Lemma 3 or Lemma 11. 

Lemma 17 If C has rotation symmetry, then there exists n G N such that 
cp n (t) = t. 

Now we can identify the parameters of ip(t) with the elements of the matrix 
A= 

U 5 J 

and one can check that the parameters in <p n {t) = (tpo oip)(t) can be 
identified with the coefficients of A n . 

Lemma 18 If C has rotation symmetry and 5 = 0, the matrix A has two real 
and distinct eigenvalues (as a consequence, it is diagonalizable over the reals). 
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Proof. When 5 = 0, the characteristic polynomial of A is A 2 — a\ — ^7. The 
roots of this polynomial are 

_ a ± V« 2 + 4/3 2 7 2 
A — - 



Since /3 • 7 7^ (because otherwise a8 — /3j = 0), the discriminant of the 
equation is strictly positive. □ 

Lemma 19 IfC has rotation symmetry and 5 = 0, then the eigenvalues Ai, A 2 
of A satisfy Ai = — A 2 . 

Proof. From Lemma 17, one may see that for some n e N, it holds that 
A n — /j, • I, with /j, e ffi. So, A n has obviously just one eigenvalue, However, 
since the eigenvalues of A n are the n-th powers of Ai,A 2 , we deduce that 
A™ = Ag. So, Ai = £ • A 2 where £ is an n-th root of the unity. But since from 
Lemma 18 both Ai, A 2 are real and distinct, we have that £ = —1. □ 

Proposition 20 If C has rotation symmetry and 5 = 0, then a = 0. As a 
result, ip(t) = k/t. 

Proof. The trace of A is a + 5. Since the trace is similarity-invariant, and 
from Lemma 19 we have that Ai = — A 2 , we get that a + 5 = 0. Since 5 = 0, 
a = 0, too. n 

Hence, in this case we have 

z(k/t) = z + e w (z(t) - zq) (20) 



and differentiating with respect to t, and solving for e ld , 

jo _ gwg • (-k/t 2 ) 

" ~ z'(t) 



(21) 



Taking modules, we get 
z'(k/t) ■ (-k/t 2 ) " 



1 = (22) 



z'(t) 

Lemma 21 Under our assumptions, the equation (22) cannot be an identity. 

Proof. If (22) is an identity then by substituting k = we get \z'(t)\ = 0, in 
which case C is not a curve. □ 

So, let £i(k) be the gcd of the coefficients of the numerator of (22) (considered 
as a polynomial in t, with coefficients in k). From Lemma 21 we know that 
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it is not identically 0. Also, let ^{k) be the gcd of the coefficients of the 
numerator of the derivative with respect to t of the right-hand side of (21). 
As in other cases, the conditions £i(k) = and ^(k) = 0, taken together, 
guarantee that (21) holds. Also, let rj(k) = gcd(£i(&;), ^(k)). Then we have 
the following result. 

Theorem 22 If C has rotation symmetry and tp{t) = k/t, then i](k) = 0. 
Conversely, if k G R fulfills r)(ko) = and no denominator in (19), (21), 
vanishes identically, then C has rotation symmetry with ip(t) = k /t, and the 
elements of this symmetry can be computed from (20) and (21). 



5.2 Case 5^0 



In this case we can assume 5 — 1. Now differentiating (19) with respect to t, 

and taking into account that (fi'it) = — -, we get 

(7* + o) 2 



z'{ V {t))-A = e ie -z'{t){ 1 t + 5) 2 (23) 

(where 5 = 1). Evaluating at t — 0, we deduce 

z\f3) ■ A = e ie ■ z'{0) (24) 



Since A ^ and because of the hypotheses, one may check that z'{0) ^ 0. 
Now by differentiating (23) again, substituting t — and taking (24) into 
account, we have that 



7 = 



2z'(0) 



z"((3)z'(0)A 

" W) 



z"(0) 



(25) 



Let us consider the real and imaginary parts of 7: 



M7) = Re (Si) ■ A - Re (|3g) 
Im( 7 ) = Im(ga).A-Im(|»), 



and let us write f((3) = Re(7), g(/3) = Im(7). Since 7 can always be taken real, 
we have 7 = f(/3), and g(/3) = 0. Now let us denote hi(/3) = Im ( Iz^p) ) ' w = 
Im ( iJfX ) ■ Observe that hi(/3), considered as a polynomial in f3, cannot be 
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identically 0. Indeed, one may see that this happens iff x'(t)y"(t) = x"(t)y'(t), 
i.e. 

x"(t) _ y"{t) 

At) At) 



for every t. Integrating twice, this implies that C is a line, which is excluded 

w 

by hypothesis. As a result, A = h((3) = . Hence, plugging this into f((3), 

hi(p) 

we get that 7 can be written as a rational function of (3, and since A = a — ^7, 
so is a, i.e. a = A + (3f(/3). Furthermore, from (24) we obtain 

A 2 = \AQ)\ 2 



So, we the deduce that: 

k'(o)l 2 



= 



(26) 



Lemma 23 Under our hypotheses, the equation (26) cannot be an identity. 
Proof. Since z' — x' + iy', then (26) can be written as 



4(x /2 + ^ 
(-x"y' + x'y") 



C 



with C a constant. If C = then C cannot be a curve. So, assume that C^O, 



and let k — ^. Then we have that 



(-x'V + x'y") 2 

( x /2 + ^2)3 



A: 



/2 



Writing tt = Jj, fci = (notice that fc > 0), we get 
k\x' 



u' 



(l+ U 2)3/2 



Integrating, we have that 
u 



hx + k 2 = 



n" 
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Now since u = — n after substituting this and taking squares, we get (after 
some calculations), 

(kix + k 2 ) 2 x' 2 _ , 2 
1 - (hx + k 2 f ~ V 

Taking square-roots, and integrating, it holds that 

^Ji-(hx + k 2 y = y + k 3 

ki v 



From here one may easily get that the implicit equation of C is that of a 
circumference, which is excluded by hypothesis. □ 

So, let be the numerator of (26), which, from the above lemma, is not 

identically 0. We also consider the equation 

d ( z\ V (t)). V '(t) \ = 
dt \ z'(t) J 



(27) 



which is obtained after solving for e 10 in (23), and differentiating with respect 



to t. Let £, 2 ((3) be the gcd of the coefficients of the numerator of — 




As in other cases, the conditions = 0, G(/3) = 0, taken together, guar- 

antee that (23) holds. Finally, let rj(/3) = gcd(^i(/3), ^ 2 (/3)), which is not iden- 
tically because of Lemma 23. 

Theorem 24 If C has rotation symmetry and (pit) = — — with 5^0, then 

7* + o 

it holds that: (i) 5 — 1, 7 = /(/?), ct = h((3) + (3 f ((3) ; (ii) i]((3) = 0. Conversely, 
ifftei fulfills: (i) i]((3 ) = 0; (ii) no denominator in (19), (25), (26), (27), 
vanishes identically, then C has rotation symmetry with the considered (p(t), 
and the symmetry elements can be obtained from (19) and (2\). 



5.3 Full Algorithm 



The following algorithm follows from the preceding subsections. 

Algorithm Rotation Symmetry: Given a curve C by means of a proper 
parametrization z(t) = x(t) + iy(t) (in complex form), well-defined at t — 0, 
such that z'(0) 7^ 0, the algorithm checks if it has rotation symmetry, and 
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computes the rotation center and the angles in the affirmative case. 



(1) ((fit) = k/t) Check if r)(k) has any real roots ki, . . . , k r G R in the condi- 
tions of Theorem 22. In the affirmative case, let z be the rotation center 
and let Q-y, . . . , 9 r be the angles they correspond to. 



conditions of Theorem 24. In the affirmative case, let zq be the rotation 
center and let Q\,...,Q S be the angles they correspond to. 
(3) If (1) and (2) have not succeeded, then return "The curve has not rotation 
symmetry". Otherwise, return the rotation center, and the list of angles. 

Lemma 21 and Lemma 23 guarantee that this algorithm terminates. The cor- 
rectness follows from the preceding subsections in this section. 



6 A Special Case 

In this section we consider a curve C parametrized by a proper parametrization 



where p(t),q(t) are polynomials, r, s G N with either r > or s > (in the 
sequel, without loss of generality we will assume that r > 0), and gcd(p(£), t 2 + 
1) = 1, gcd(q(t),t 2 + 1) = 1. This subclass of rational curves is important 
because it contains most of the so-called trigonometric curves. These are the 
real plane curves where each coordinate can be given parametrically by a 
trigonometric polynomial, that is, a truncated Fourier series: 



with aj,fej,Q,dj G R. These curves appear in various areas of mathematics, 
physics and engineering (one may see the introduction to [5] for a list of dif- 
ferent contexts where they arise). They have been considered for example in 
[1], [5], [12], [23], [24], [13]. From the results in [5] (see Theorem 2.1 therein), 
it follows that every trigonometric curve has either a simple trigonometric 
representation (i.e. a parametrization of the above kind, which is injective for 
almost every point when restricted to 6 G [0, 2ir]), or a polynomial simplifica- 
tion, i.e. a polynomial parametrization which coincides with the curve when 




Check if rj(/3) has any real roots (3i, . . . , (3 S in the 




x = Er=oK cos (^) + b k sm(k9)) 
V = J2'k=o( c kcos(k8) + d k sin(k9)) 



18 



the parameter is restricted to a certain real interval [a, b}. For the first ones, 
De Moivre's formula together with the usual parametrization 

cos(0) = sm(V) = 



1+t 2 ' v ' 1 + t 2 



applied on the simple representation provides a parametrization of the kind 
considered in this section. 

In the sequel, we will denote m = deg(p(t)), n = deg(q(t)). Now for these 
curves, and specially for central symmetry and mirror symmetry, we can pro- 
vide a sharper class of temptative <^(t)'s, which shortens (in some cases, dra- 
matically, as it is shown in the next section) the computations. 



6.1 Central Symmetry 



Following the ideas in Section 3, if C has central symmetry then z(tp(t)) = 2zq — 
z(t). If we consider only the first coordinate, we get that x(<p(t)) = 2x — x(t); 
substituting ip(t) and performing some calculations, we have that 

m = 2x _ Pit) 



where P(t) is a polynomial which can be written as P(t) = a m {at + (5) m + { r yt+ 
S)P(t). Notice that since gcd(p(t),t 2 + 1) = 1 and t 2 + 1 is irreducible then 
(at + (3) 2 + (yt + 5) 2 cannot divide P(t). Additionally, from the irreducibility 
of t 2 + 1 we deduce that: 

(a) If m > 2r, then yt + 8 must divide t 2 + 1. Hence, the only possibility is 
7 = 0. Furthermore, (at + (3f + (7* + 5) 2 = {7 = 0} = a 2 t 2 + 2a/3t + /3 2 + 5 2 
must also divide t 2 + 1. So, a/3 = and a 2 = /3 2 + 5 2 . Since (otherwise 
we would get a5 — ^7 = 0) we get /3 — 0, a — ±5. Discarding the solution 
a = 5 (which implies x(t),y(t) constant) we get a = —5, i.e. (p(t) = —t. 

(b) If m < 2r, then (at + /3) 2 + (-ft + 5) 2 = (a 2 +-f 2 )t 2 + 2(a/3+^5)t + l3 2 + 5 2 must 
divide t 2 + 1. So, a (3 + 7<5 = and a 2 + 7 2 = (3 2 + 5 2 . Now we distinguish 
two cases: 

■ If 5 = 0, then a/3 — 0. If /3 — then we get a 2 + 7 2 = 0, and since a, 7 are 
real, a = 7 = 0, which contradicts a5 — ^7 7^ 0. So, in this case a — 0, 
7 = ±/3, and hence tp(t) = 1/t or (p(t) = —1/t. 

■ If 8 ^ 0, then we can assume 8 — 1. Moreover, from Subsection 3.2 we 

know that a = —5, and so a = —1. Hence, 7 = ft and (p(t) = ~*~ — . 
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As a consequence, the only possibilities for ip(t) are: —t, 1/t, — 1/t and ip(t) = 
-t + /3 

Pt + 1" 



6.2 Mirror Symmetry 



The same idea can be used in the case of mirror symmetry and rotation sym- 
metry. In the first case, the condition (5) can be written as 

zjtj = z^+ [z(<p(t)) -z ] -e 2ie 



Decomposing z(t) into real and imaginary parts, we get 

p(t) t R(t) ■ cos(2#) S(t) -sin(20) 

— x n + 7 ; ; ?m ft ; ; tw^, ; 7 ; ; r\m„ + 



(t 2 + l) r u (t* + 5) m - 2r ■ [(at + (3) 2 + (t* + 5) 2 } r ( 7 t + 5) n ~ 2s ■ [(at + (3) 2 + ( 7 i + 5) 2 Y 
q(t) _ t M (t) • sin(20) N(t) ■ cos(20) 



(t 2 + 1)* yu ( 7 t + 5) m - 2r • [(at + /3) 2 + ( 7 t + 5) 2 } r ( 7 t + 5)™- 2s • [(at + /3) 2 + ( 7 t + 5) 2 ] s 

where -R(t), S(t), M(t), N(t) are polynomials, = x — a; o cos(20) + y sin(26), 
y$ = —y — y cos(29) — x sin(29). Now arguing as in Subsection 6.1, we have 
that: 

(1) If m — 2r < and n — 2s < simultaneously, then the only possibilities are 
<p(t) = 1/t, <p(t) = -1/t, or <p(t) = 

(2) Otherwise there is one more possibility, namely ip(t) = —t. 

So, the possible </?(i)'s are exactly the same as in the case before. 



6.3 Rotation Symmetry 



Arguing as in the above subsection, if m — 2r < and n — 2s < then in 
the case 5 = 0, the only possibilities are ip(t) = 1/t or ip(t) = —1/t. However, 
when 5 ^ 0, since in this case it is not true in general that a = —5, we cannot 
reach a simple form for <p(t), similar to that in the preceding subsections. 



7 Implementation and experimentation 



From the algorithms in Sections 3, 4, 5, 6 we have that the detection of central, 
mirror or rotation symmetry reduces to checking the existence of real roots 
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of certain univariate polynomials. As a consequence, it can be performed de- 
terministically. However, in order to find the elements of the symmetry (i.e. 
symmetry centers, symmetry axes, etc.), in general a numeric approach of the 
real roots of the polynomials is required. In this sense, our algorithms have 
been implemented and tested in the computer algebra system Maple 15. As for 
the computation of the elements of the symmetry, we have observed that he 
precision needed in this computation must be increased accordingly with the 
size of the input; in particular, as the degree is increased, the infinity norm of 
the polynomials arising in intermediate steps of the algorithm grows as well. 
So, in our implementation we have followed the following strategy: 

(i) The real roots of the polynomials are approached numerically with a certain 
precision which is initially computed based on the size of some intermediate 
polynomials, appearing when certain steps of the algorithm are executed, 
which we have heuristically detected as "good indicators" . 

(ii) The symmetry elements are approached using this precision. 

(iii) Afterwards, these elements are tested to check their accuracy. For this pur- 
pose, we compute an index, which must be very close to 0: if we identify 
that this index is not small, the precision is increased, and the computa- 
tions start again. The index is computed as follows: we generate a number 
of points (typically, 100) in the curve, we compute the symmetric points of 
these with respect to the elements that we have determined, and we "esti- 
mate" the average distance from these points to the curve. In general, we 
observe a very well performance of the indicators pointed out in (i). We 
must say that as the size of the output grows, the time consumed by this 
last part (the test) is bigger and tends to be even longer than the time 
consumed in the computation of the symmetry elements; the reason behind 
is the necessity, in these cases, of a big (sometimes, huge) number of digits 
to guarantee an accurate computation. 

Also, in our experiments we have seen that generally we can manage curves up 
to degree 10 in a fast way. In some cases we can go a bit further, although this 
seems to be more difficult in the case of mirror symmetry. We can always go 
further if we avoid the computation of certain polynomials, but this implies 
that the answer yes/no on the existence of the considered symmetry is not 
deterministic anymore. Still, in this case we can always determine the temp- 
tative elements of the symmetry, and check a fortiori whether they are fake 
or not (see later the timings for the mirror symmetry algorithm) by using the 
test above. Finally, let us mention also that the timings improve dramatically 
in the case of the curves addressed in Section 6. All the experiments have 
been performed on an Intel Core revving up to 2 GHz., with 8 Gb. of RAM 
memory. 

So, let us address in more detail each of the algorithms. The following table 
shows some data corresponding to the algorithm for checking central symme- 
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try. Here, we provide the timing (in seconds), which includes the computation 
of the symmetry elements and the test on their accuracy, as well as data about 
the input curves (maximum infinity norm of the numerators and denomina- 
tors, degree of the parametrization). 



Curve 


deg t (<£(*)) 


Norm 


Sym. (Y/N) 


Timing 


Comments 


Trisectrix 


3 


3 


N 


2.372 




Lemniscata 


4 


30 


Y 


2.106 




Three-leaved rose 


4 


3 


N 


2.403 




1 


5 


9 


Y 


2.262 




Astroid 


6 


86400 


Y 


2.918 


71 digits 


2 


7 


108 


Y 


3.775 


58 digits 


3 


9 


800 


Y 


11.482 


219 digits 


4 


11 


460 


Y 


37.128 


276 digits 


5 


14 


22801 


Y 


692.458 


638 digits 


16-leaved Rose 


18 


25740 


Y 


2.465 





For degrees higher than 10 we get worse timings, in general, due mainly to 
two causes: (1) the cost of computing the function /(/?) (i.e. the relationship 
7 = /(/?)); (2) the cost of performing the test (because of the required number 
of digits). The curve 5 illustrates this situation. Still, in certain cases and 
even in presence of high degrees, we may get a good performance when a 
symmetry center is detected without computing f((3) (the case of the 16- 
leaved rose). Furthermore, detecting curves with no central symmetry is faster. 
Now the following table shows data corresponding to the algorithm for mirror 
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symmetry. 



Curve 




Norm 


Sym. (Y/N) 


Timing 


Comments 


Descartes' Folium 


3 


3 


Y 


3.541 




Epitrochoid 


4 


288 


Y 


3.042 




Lissajous Curve 


8 


20 


Y 


3.573 




Offset Cubic Curve 


8 


54 


Y 


27.097 




Offset Cardioid 


8 


78732 


Y 


17.176 




9 


10 


99 


Y 


35.085 


61 digits, (*) 


fO 


12 


6 


Y 


8.674 


w 


ff 


13 


371 


N 


40.046 


19 digits, (*) 


12 


16 


9 


Y 


192.194 


17 digits, (*) 


Sixteen-leaved Rose 


18 


25740 


Y 


65.941 


28 digits, (*) 


Twenty-leaved Rose 


22 


369512 


Y 


168.824 


34 digits, (*) 



In our experiments, we have observed that most of the time is spent in the 
computation of the expression (18), and in the test of the symmetry axes. In 
fact, the timing for the last six curves (those with a (*)) corresponds to the 
computation without explicitly determining (18). This implies that the set of 
symmetry axes to be tested at the end of the algorithm may include some 
which is fake; however, in our experiments these fake axes have been always 
identified by means of high values in the index (which is computed in the 
test). Now the following table corresponds to the algorithm to detect rotation 
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symmetry. 



Curve 


deg^(t)) 


Norm 


Sym. (Y/N) 


Timing 


Comments 


Trisectrix 


3 


3 


Y 


4.399 




Epitrochoid 


4 


288 


N 


4.352 




Three-leaved-rose 


4 


3 


Y 


4.992 


36 digits 


Astroid 


6 


8 


Y 


11.123 


23 digits 


Four- leaved-rose 


6 


12 


Y 


6.209 




Lissajous Curve 


8 


20 


N 


4.399 




13 


9 


5 


N 


178.466 


23 digits 


Eight-leaved rose 


10 


140 


Y 


26.957 


28 digits 


Twelve-leaved rose 


12 


1848 


Y 


129.246 


42 digits 


Sixteen-leaved rose 


16 


25740 


Y 


510.701 


55 digits 



Finally, in the following table we compare, for some curves defined by a 
parametrization of the kind studied in Section 6, the timings correspond- 
ing to the algorithms for checking central and mirror symmetry in Sections 
3 and 4 (denoted "ord"), and the implementation of the ideas in Section 6 
(denoted "spec"). One may notice that the latter outperforms the former by 
a considerable factor. 



Curve 




Norm 


Sym. tested 


Timing (ord.) 


Timing (spec.) 


17 


25 


504 


C 


> 1000 


1.966 


14 


20 


504 


C 


> 1000 


1.482 


15 


26 


3696 


C 


> 1000 


2.606 


Sixteen-leaved rose 


18 


25740 


M 


139.886 


33.322 


Twenty-leaved rose 


22 


369512 


M 


502.074 


53.867 


14 


20 


504 


M 


> 1000 


2.543 


17 


25 


504 


M 


> 1000 


1.654 
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8 Appendix I: Proof of Lemma 14 in Section 4.2 



In order to prove Lemma 14 we need several previous considerations and 

z"(0) 

results. First let a,b G R fulfill a + ib = . Taking into account that 

z'(0) 

z = x + iy, one may see that 

_x'(0)x"(0)+y'(0)y"(0) 
a — 



b = 



x' 2 (0) + y' 2 (0) 
x"(0)y'(0)-x'(0)y"(0) 

x /2 (0) + y' 2 (0) 



Furthermore, from the expression (16) in Section 4.2 we have that (here we 
use the above notation): 

(2 + at){x'x" + y'y") + bt(x'y" - x"y') + 2a(x' 2 + y' 2 ) + t(x" 2 + y" 2 ) 

~ (tx" + 2a;') 2 + (ty" + 2y') 2 

(2 - at)(x'y" - x"y') + btjx'x" + y'y") + 2b(x' 2 + y' 2 ) 

~ (tx" + 2x') 2 + {ty" + 2y') 2 

For simplicity, in the above expressions we have written t where we should 
have written (3; we will proceed in the same way along the appendix. Now 
notice that since 7 can always be taken real, from the above expression for 
1111(7) we deduce that 

(2 - at)(x'y" - x"y') + bt(x'x" + y'y") + 2b(x' 2 + y' 2 ) = (28) 



is an identity. 

Lemma 25 Under our hypotheses, b 7^ 0. 

Proof. If b = 0, then substituting in (28) we get x'y" — x"y' = 0, which leads 
to y' = cx', i.e. y = cx + d; so, C is a line, which is excluded by hypothesis. □ 

Lemma 26 There does not exist t e C such that x" (to)y' (to) — x' (to)y" (to) = 
0. 

Proof. If such &to £ C would exist, we might reparametrize the curve with u = 
t+t . However, the resulting parametrization would fulfill 6 = 0, contradicting 
Lemma 25. □ 

Then the following result on the choosing of a, b holds. 

Lemma 27 Without loss of generality we can assume a = 0, 6 = 2. 
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Proof. We proceed in a constructive way. We want to find a reparametrization 
of the form 

. . m ■ u + n 

t(u) = — 

v ' r-u + 1 

with m,n,rei,m-nr^O such that x = x(t(u)),y = y(t(u)) fulfills that 
a = 0, b = 2. So, let 

x'x" + y'y" 
a = a (i ,y){u) = _ /2 - — 2 



One can check that 

x'x" + y'y" _ x'{t)x"{t)+y'{t)y"{t) , f_ 
x' 2 + y' 2 ~ x' 2 {t) +y' 2 {t) 1 + t 7 



Furthermore, i! = ™ u ™ )2 and t" = -2r ( ™ + "^ . Substituting -u = 0, we get 

x'(n)x"(n) +y'(n)y"(n) / . n 

a = ,2) L a ( \ > {m-nr)-2r 29 



Similarly, 

~ y>(t)x"(t)-x>(t)y"(t) 



and substituting in u = 0, we have 

6 = ^?!'y (n) (m - nr) (30) 



So, we want to find m,n,r e K such that a = 0, 6 = 2. By imposing these 
conditions and dividing (29) and (30), we obtain 

x'(n)x"(n) + y'(n)y"(n) 

T = 

y'(n)x"(n) — x'{n)y"{n) 



Notice that by Lemma 26, the denominator of the above expression does not 
vanish for any real value of n. Furthermore, from this equality and (30) it 
holds that 

2(x'(n) 2 + y'(nf) 

m — nr = — — — — — — — 

y'(n)x"(n) — x (n)y "(n) 
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So, by choosing no G K such that the above expressions for r and m — nr are 
well defined, and m — n r 7^ 0, a reparametrization with the desired properties 
is found. □ 



Now we can prove the following result. Recall here the notation 7 = /(/?) 
(however, we will write f(t), instead). 

Lemma 28 IfC exhibits mirror symmetry and £1,^2 are identically zero, then 
for every t EC, it holds that \z'(t)\ 7^ 0. 

Proof. If £2 is identically 0, then from the expression (17) we have \z'(t)\ = 
u(t), where u(t) is a rational function. Writing z'(t) = x'(t) + iy'(t), there 
exists an analytic function v(t) fulfilling x' — u • cosf , y' = u • sint>. Using the 
expression given for 7 in the beginning of the appendix, and taking Lemma 
27 into account, one can see that 

_ 2uu' + 2tu 2 v' + t((u / ) 2 + u 2 (v / ) 2 ) 
/U Atuu' + 4u 2 + t 2 ((u') 2 + u 2 (v>) 2 ) 



Additionally, since Im(7) = 0, from the expression for 1111(7) a ^ the beginning 
of this appendix we get that u 2 v 1 + tuu' + 2u 2 = 0. By combining this with 
the above expression for fit), we have that 

= ~2u + tu' 



Now let us argue by contradiction. For this purpose, let to G C be a zero of 
|z'(t)|. If t — 0, then from the expression (17) in Subsection 4.2 we deduce 
that \z'(t)\ is identically 0, and hence C is not a curve. So, t 7^ 0. Also from 
(17) one may see that if t = t is a zero of |-z'(t)|, then it must be a pole of f(t), 
which is a rational function. Now since t = to is a zero of u(t) (of multiplicity 
m), then u(t) = (t — to) m r(t), where r(t) is rational and r(to) 7^ 0. Hence, 
differentiating this expression for u(t), substituting in the above expression 
for f(t), and eliminating common factors in the numerator and denominator, 
we get that 

u' _ mr(t) + (t - to)r'(t) _ 1 

im '^'° ~ 2u + tu' ~ im ^*° ~ 2(t - t )r(t) + mtr(t) + t(t - t )r'(t) ~ " ' " ~J~ 



Since to 7^ 0, we see that the above limit is not equal to infinity, and therefore 
t — to is not a pole of fit), which is a contradiction. n 

In the sequel, we will use the notation u(t) = \z'(t)\, x' — u • cost" , y' — u- smv 
introduced in the proof of the above lemma. Now we are ready to prove Lemma 
14. 
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Proof of Lemma 14. If v is constant then C is a line, and we have finished. 
Also, if u is identically then C cannot be a curve. So, in the sequel we assume 
that v is not constant, and u is not identically 0. From Lemma 27, it follows 
that £i being identically is equivalent to the following identity: 

2(x'y" - x"y') + 2t(x'x" + y'y") + 4(x' 2 + y' 2 ) = 



With the above variables u, v, this equation is transformed into 

v' = -t 2, 

u 

which in turn yields 

J V -dt = -m(C* • t 2 u) 

Let R(t) = Cq ■ t 2 u. Since R(t) is rational, we can write 

R(t) = C -(t- to)" •••(<- t p ) n " ■ {t 2 + a t + b ) mo ■ ■ ■ (h + a q t + b q ) m « 

with rii G Z, rrij G Z, C G K and a 2 — 4^- < 0, Vj G {0, 1, . . . , g}; furthermore, 
all the tj's are different, and all the terms t 2 + ajt + 6/s are also different. 
So, differentiating the logarithm of the above expression and doing some easy 
computations, we have that 

n t : : n p tp a t + 2b a g t + 2b q 



v' = n H h n p + 2(m H h m q ) + - — — H h —^7 — m o 



t-t ' ' f-*„ " H) t 2 + a t + b mq t 2 + a q t + b 



Since we are assuming that v is not constant, we get that 

y' 

v = arc tg— = I + h + h + h 
x' 



where 



Iq — C*2 + Cit 

h = \n\t - t \ not ° ■ ■ ■ \t - tpl™^ 
h = _zMooi n |i2 + flot + 6() | »ln|t 2 + a,t + 6, 

^a r c tg (^-...--#^- tg 
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Now observe that v and I3 are necessarily bounded as t moves in K. So, noto = 
■ ■ ■ — n p t p = 0. Since all the t^s are different, we have that at most one of 
them can be equal to 0, and therefore there must be at least p elements in 
{no, ni, . . . , Up) which are 0. W.l.o.g., let us assume that n\ — . . . — n p — 
0. For the same reason (i.e. v and I 3 bounded for t G K), C\ = 0. As a 
consequence, we deduce that 

, . n( x a t + 2b a q t + 2b q 

v = n + 2(m H h m q ) - m — — — m i7Ts TVT 

t 2 + a t + b t 2 + a q t + b q 

where no + 2(m + • • • + m q ) — C\ — 0. Substituting t = we get that 
v'(0) = n . However, since v' = —t— — 2 we also deduce that v'(0) = —2 
(notice that u(0) 7^ because we assumed ^'(0) 7^ 0, see the beginning of 

Section 4), and therefore n = —2. As a result, we get that m H Vm q = 1, 

and 

. a t + 2b a Q t + 2b q 

v = —m — • • • — m q 



t 2 + a t + b q t 2 + a q t + b q 

This means that 

j = (in (C ■ t~ 2 ■ (t 2 + a t + b ) mo ■ ■ ■ {t 2 + a q t + 6,)^))' 

Therefore, 

u = C (t 2 + aot + b )- m » ■■■(t 2 + a q t + b q )- m " 

Now since all the m^s are integers, and m + • • • + m q = 1, there are two 
possibilities: (a) all the m^s are except, say, m = 1; (b) some m/s are 
positive, and others are negative. However, (b) cannot happen because this 
would contradict Lemma 28. Hence, u = Co(t 2 + aot + bo)^ 1 , and 

, a t + 2b 

v = —m - 



t 2 + a t + b 

Nevertheless, do = because otherwise the primitive of v would contain a 
logarithm. So, finally we get u = — t 21 ^ ho C and v' = — t 2 ^ . So, u = kv'. 
Hence, we get 

x' = kv'cos(v) 
y' = kv'sm(v) 

Integrating, we deduce that C is a circle. □ 
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